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SUMMARY 

The theory, design, and effect on cooling of the constant- 
velocity axial-flow fan as related to the problem of high-altitude 
aircraft engine cooling are discussed in the present report. The 
existing theory of the axial-flow fan, including a discussion of 
the fan losses, is summarized. Nondimensional design and perform- 
ance formulas are developed for stator-rotor , rotor-stator, and 
rotor-alone fan combinations. Charts based on these formulas are 
given to facilitate performance estimates and design. Some illus- 
trative fan designs are made for an altitude of 40,000 feet. 
Finally, the effect of a fan on cooling power is quantitatively 
discussed and illustrated. 

The results of the present study indicate that an axial-flow 
fan may be used to give a moderate pressure rise to a large flow 
quantity of air; that existing theory provides a satisfactory basis 
for the design of axial-flow fans; that the use of stationary guide 
vanes and small rotor tip clearances is essential for efficient fan 
design; that further research is needed for the development of air- 
foil sections especially suitable for use in high-speed, high- 
solidity rotors; that the losses in internal cooling power at high 
altitude are large when the cooling is marginal; and that an axial- 
flow fan more than pays for itself in cooling power saved when 
cooling conditions are marginal.. 



INTRODUCTION 

As part of a general study of various aspects of the problem 
of high-altitude engine cooling (see reference 1), the theory, 
design, and effect on cooling of the constant-velocity axial-flow 
fan are considered in the present report, 
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An axial-flow fan may- be effectively used as an engine -cooling 
aid at high altitude, where the necessary cooling pressure is not 
entirely available or where, if obtainable with large exit flaps 
and flap deflections, the resulting drag would be prohibitive. 
Although the fan itself acts primarily to increase the pressure 
of the air passing through it, several other beneficial effects 
result from its use in an airplane. 

First, the pressure increase of the cooling air allows a decrease 
in the ultimate kinetic energy of the cooling air following the air- 
plane. The net power saved in this way may be more than the power 
required to drive the fan, with a resulting increase in propulsive 
efficiency for the airplane as a whole* This effect will be dis- 
cussed quantitatively in a later section. 

Furthermore, an axial-flow fan mounted at the entrance of a 
cowling can decrease the lower limit of cowling inlet velocity 
for smooth-entry flow and improve the axial velocity distribution 
around the inside of the cowling. The reason probably lies in the 
fact that the pressure rise across a fan goes up if the flow through 
it goes down and vice versa; hence, any variations in the axial 
velocity distribution in front of the fan produced by boundary 
layers or by the attitude of the airplane tend to be ironed out by 
suitable pressure gradients produced by the fan. It may also be 
mentioned that a fan at the entrance of a cowling is subject to 
the same angle -of -attack effect 3 as a propeller; that is, if the 
fan axis is tilted upward relative to the incoming flow, the blades 
on the downgoing side are more highly loaded, and hence produce a 
higher pressure rise, than the blades on the upgoing side. 

The design of an axial-flow fan must be carried out with the 
characteristics of its power supply in mind. There are three pri- 
mary sources of power on an airplane in flight; namely, the engine, 
the engine exhaust gases, and the motion of the airplane relative 
to the air. The corresponding devices for transferring this power 
to the fan would be a gear or direct drive from the engine, an 
exhaust-gas turbine, and a windmill. Each device can be used 
practicably to drive the fan and the choice therefore depends on 
considerations of engine and cowling design. 

The theory of the constant-velocity axial-flow fan operating 
in an incompressible fluid, together with design formulas and 
charts, is presented in the sections immediately following. The 
design procedure as well as the cooling possibilities at altitude 
of the axial-flow fan is then illustrated by the design of three 
arrangements of single-stage axial-flow fan for operation at 
40,000 feet. The concluding section deals with the internal - 
power losses in a cooling installation, with compressibility 
effects included. 
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SYMBOLS 

C[ dynamic pressure {^^^j 

p static pressure 

Ap change of static pressure 

V velocity 

V free- stream velocity 
p density 

G mass flow per unit area (pV) 

Q quantity rate of flow 

M mass rate of flow 

H heat rate of flow 

P power 

r\ efficiency 

s cascade spacing 

w relative velocity 

X x-directicn (tangential) force per blade per unit span 

V y -direction (axial) force per "blade per unit span 
Z "lift" per "blade per unit span 

d drag per blade per unit span 

T circulation 

c blade chord 

Cj section lift coefficient 

section drag coefficient 

Cp section pressure coefficient; also, specific heat at 
constant pressure 
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c v specific heat at constant volume 

c T section torque coefficient 

r radial location of blade element 

x nondimensional blade-element radial location (r/R 0 ) 

E Q fan radius at outer section 

R^ fan radius at hub section 

00 angular velocity 

u velocity in rotational direction 

n rotational speed, revolutions per unit time 

B number of blades 

N rotation parameter (tor/Vf) 

7 circulation parameter (p/cVf); also, ratio of specific 
heats (Cp/c v ) 

G solidity (Bc/2*r) 

T torque; also, temperature 

T v average cooler temperature 

a angle of attack 

p blade angle measured from zero lift line to fan plane 

(3 1 blade angle measured from chord line to fan plane 

0 relative stream angle measured from relative velocity direction 
to fan plane 

f Collar cascade lift factor 

6 Ruden cascade thickness correction 

5 blade profile area 

D diameter; also, reciprocal of profile aspect ratio (S/c 2 ) 
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a o 


slope of section lift curve 


a *o 


angle of attack for zero lift 


A 


cross-sectional area 


a 


speed of sound 


Subscripts: 
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p 


pressure 


X 


x direction 


y 


y direction 


1 


entrance 


fi 


upstream of fan 


tz 


dovnstream of fan 


z 


cooling-resistance entrance 


3 


downstream of cooling resistance 


3e 


cooling-resistance exit 


4 


ultimate exit condition 
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s 



stator; also, stagnation 



t 



tip clearance 



max 



maximum 



THEORY OF TEE AXIAL -FLOW FAN 



The axial -flow fan consists essentially of one or more sets of 
rotating impeller vanes , hereinafter called "rotors," which may "be 
operated in conjunction with one or more sets of stationary guide 
vanes, hereinafter called "stators." (See fig. 1.) A rotor operating 
alone increases the static pressure and puts kinetic energy of rota- 
tion into the fluid passing through it in the axial direction. A 
stator may "be used behind or downstream of the rotor to convert the 
kinetic energy of rotation into additional static-pressure increase • 
There is then a static-pressure ri3e across both rotor and stator. 
A stator may, on the other hand, be placed ahead or upstream of the 
rotor as in figure 1. In this case the stator twists and reduces 
the static pressure of the flow, which is then straightened and 
increased in pressure by the rotor. 

The analysis of the axial-flow fan can be greatly simplified 
if three effects are neglected. The first is concerned with the 
various boundary- layer motions, the so-called secondary flows, 
which arise in response to the pressure gradients set up by the 
fan operation. So far only a qualitative insight into these 
phenomena has been gained. (See reference 2, sec. V.) The second 
effect is that of the radial velocities produced by nonunif ormity 
of pressure increase over the fan area. Ruden (reference 2) has 
investigated this effect and gives a method of averaging the fan 
blade-element characteristics calculated for zero radial velocities 
for the determination of the over-all fan characteristics. At the 
design operating point the design criterion is uniform total pres- 
sure, with the result that radial velocity components and a nonuni- 
form energy distribution in the flow through the fan do not arise. 
Eadial velocities may therefore be neglected if the design operating 
conditions are adnered to. For other operating conditions, at which 
radial velocities may not be negligible, the averaging method of 
reference 2 can be used. Finally, the spanwise mutual interference 
of the blade elements is neglected. This assumption is well justi- 
fied when the blade circulation is constant along the span (the 
design condition) or when the number of blades is large, as in an 
axial fan designed for high pressure rise. 
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Consider the annular area of a set of "blades (fig. 2) included 
"between the radial stations r and r + dr and rotating with 
angular velocity CO in a flow that has "both an axial velocity 
component Vf and an initial rotational velocity component u^. 
Positive directions are indicated "by the arrows in the figure. 
Under the assumptions made, the annular set of "blades may "be 
regarded as a two-dimensional cascade of airfoils. It is desired 
to express the forces on an airfoil of the cascade in terms of a 
lift and drag force analogous, if possible, to the forces on an 
isolated airfoil. To this end the cascade is regarded as station- 
ary, its original velocity being imparted in reverse to the 
fluid. The equation of continuity, Bernoulli's theorem, and the 
momentum theorem are then applied to the fluid within the boundary 
ABCD3F where ABC and FED are relative streamlines and AT 
and CD are sections X and 2 parallel to the cascade direction. 
Drag effects are for the moment neglected, with the result that 
the flow is uniform at the sections 1 and 2. 

The equation of continuity is simply 

V f = Constant (l) 

Bernoulli ! s theorem for sections 1 and 2 is 



PjL + |p j\r f 2 + (cor 4- u^ 2 = p 2 + |p j\ f 2 + (tor + u 2 ) 2 j 



(2) 



The' momentum theorem gives the force components X and Y on the 
blade element per unit length of blade span. If it is remembered 
that the pressures at similar, points on the streamlines ABC and 
FED are the same and if the continuity condition of equation (l) 
is used, then the momentum equations are 



x direction (tangential) 
X = psV f >(u- L - Ug) 

y direction (axial) 
Y = s(p 1 - p 2 ) 
or, using Bernoulli's equation (2), 



(3) 



Y - ps! u 2 



\ / Ut + Ug^ 



( V 

^ u 2 " *l) ^ 



tor + 



J 



(4) 
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The forces X and Y are the component a of a "lift 11 per unit span 

I = pvr (5) 

where 

r = s(u 1 - u 2 ) (6) 

and 



v? - 



The cascade spacing s of equations (3) and (4) has "been expressed 
in equation (7) in terms of the circumference and number of "blades 

2*r 
3 = ~ 

The quantity F defined by equation (6) is the circulation 
around the blade element. The velocity v is the mean of the 
relative velocities upstream and downstream of the annulus. For 
large values of the ratio of blade spacing to blade chord s/c, the 
velocity w approximates the actual flow in which the blade element 
acts. The lift I on the blade element acts at right angles to the 
mean relative velocity w, in analogy with the lift on an isolated 
airfoil. In addition to the lift there is a profile drag d, which 
is customarily assumed to act in the direction of w whereas its 
magnitude must be given by experiment. Although this assumption 
is in error , as discussed in the section "Fan Losses/' it will be 
at least temporarily retained because of its simplicity and because 
the over-all error is not large. 

With the aid of equations (5) and (2), the lift, the drag, and 
the mean relative velocity can be expressed nondimensionally . If 
c is the local blade chord, 



c 



1_ = 42L (8) 



1 J pw 2 c ~ v/V f 



Za - 



(9) 
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'f 



where 

the rotation parameter 
the circulation parameter 
the solidity 



N = 



7 = 



or 



r 

cV, 



Be 

and the mean relative stream angle <f) is shown in figure 2. The 
pressure rise Ap and the torque T per unit length spanwise of 
the rotor are obtained by resolving these forces parallel and 
perpendicular to the fan axis. The nondimensional results are 



Ap f 



Crp 



2*r 2 q f 



=^ + ^t( B+ ^"^) (12) 



The power output dP Q of the annulus is. dQAp where 

dQ » 2jrrV f dr 



Thus , 



dP Q = 2*rV f Ap f dr (13) 
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and the pover input to the annulus is 



(14) 



which gives a section efficiency of 




(15) 



These quantities, as veil as the pressure and torque, may be suitably 
integrated over the fan area; but due account should be taken of the 
tip-clearance losses discussed later. 

By using equations (8) and (10) the circulation parameter y 
can be eliminated from equations (11) and (12), giving the pressure 
and torque in terms of known quantities and the lift and drag coef- 
ficients of the blade element. The problem thus resolves itself 
into the determination of the angle of attack a of the blade 
element relative to the mean flow w to give a lift coefficient 
Cj and the determination of the corresponding drag coefficient c^. 
Before these questions are taken up, the fan losses will be discussed 
and the performance formulas of various combinations of rotor and 
stator derived. 



Three main energy losses and one potential source of energy 
loss arise in the action of a rotor: (l) the profile-drag losses 
on the blades, (2) the tip-clearance losses, (3) the inner-wall 
losses, and (4) the rotational kinetic energy of the flow downstream 
of the rotor. 

(1) The profile drag of the blade element has been approxi- 
mately taken into account by a drag at right angles to the lift. 
A more detailed analysis of a cascade, similar to that previously 
given for the frictionless case but including the effect of a wake 
that starts from the trailing edge of each airfoil and passes down- 
stream with gradually increasing thickness, shows that equations (1) 
and (2) must be altered. If the subscript 2 denotes the ultimate 
downstream condition in which complete mixing of the wake has 
occurred, then equation (2) for the pressure rise has an additional, 
essentially positive term on the right-hand side. This term is a 
function of the velocities and dimensions of the wake at an arbi- 
trary downstream location, and equation (2) shows that it is invari- 
ant with respect to distance downstream from the cascade. It 
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represents a drag in the axial direction. The lift, acting perpen- 
dicular to the mean flow, is still given "by equations (5) to (7) 
but, as noted, with the subscript 2 referring to ultimate down- 
stream conditions. 

The resultant force on an airfoil in a cascade may of course 
be resolved, although somewhat artificially in view of the pre- 
ceding remarks, into components in any direction. It was found 
in reference 3 that, for a typical rotor blade, the profile-drag 
coefficient in the mean-flow direction was greater than the iso- 
lated section value by about 30 percent, and this fact will be 
used in the section "Illustrative Fan Designs." It seems, however, 
that drag measurements on airfoils in cascade should be presented 
in terms of a more rigorous analysis such as just indicated. 

(2) The tip- clearance losses are usually much greater than 
the profile-drag losses. They are caused by flow leakage from 
the high-pressure to the low-pressure side around the blade tip. 
The loss in efficiency due to tip clearance is evidenced as a 
reduction in pressure rise of the fan, the power taken by the 
fan remaining about the same for all clearances. Thi3 effect 

is probably largely caused by turbulent mixing in the axial 
direction. 

In reference 2, the tip-clearance losses were determined 
experimentally for an eight-blade rotor-stator combination with 
rotor solidity a varying from 0.79 at the hub to 0.52 at the 
tip section. The hub diameter was 9.85 inches and the outer-wall 
diameter was 19,7 inches. Figure 3 presents the mean of curves 
determined for three different blade loadings. The separate 
curves depart from the mean by about 1 percent. The reduction in 
pressure due to tip clearance is large, about 7.5 percent for a 
clearance of only 2 percent of the blade span. 

(3) The inner-wall losses are caused by the boundary layer 
on the inner wall of the rotor. By the same argument as pre- 
viously given these losses should also manifest themselves as a 
reduction in pressure rise. The dead air at the hub, however, 
is swept to the tip by the centrifugal force of rotor rotation, 
with the result that these losses may be largely included in the 
tip-clearance losses. Some measurements of wall losses in a 
stationary cascade are given in reference 4. The results are 
not directly applicable to a rotor, however, because of the span- 
wise flow of dead air just mentioned. 
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(4) The rotational kinetic energy of the downstream flov repre- 
sents a loss if it is not transformed to pressure energy "by a set of 
stationary guide vanes or by other means. This loss represents a 
drop in pressure rise and is included in equation (11) for c . 



PERFORMANCE FORMULAS FOR VARIOUS FAN COMBINATIONS 

Stator Upstream of Rotor 

The arrangement in vhich the. stator is upstream of the rotor 
will' be referred to as a "stator-rotor combination. " Subscript s 
will designate stator quantities; r, rotor quantities. 

The initial flov into the stator is assumed to have no rota- 
tional component, that is, \x± = 0 for the stator. The flov 
emerging from the stator will have a rotational velocity component 
7 a a g Vf that is to cancel the rotational velocity component 7 r O r V^ 
of the rotor. The product yo (no subscripts are necessary) is 
therefore the same for rotor and stator at any radial location; 
hence, for the stator the relative velocity v g /Vf is, from equa- 
tion (10), 



and the lift coefficient is given by 



O q c ? = — 4— 
3 V V f 



The pressure rise, vhich is negative, is, from equation (11), 

APf s 



c 



Ps 



If 



= - (7C) 2 - 0 8 C ds ^ (17) 

ttn . 

For the rotor, if it is remembered that — = yo and if the small 

effect of the stator drag on the momentum change through the stator 
is neglected, there is obtained from equations (8) and (10) to (12) 
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1 

sin 0 r 



17 



Cm = 



+ (7a) 2 - a r c d 

T 



r V f 



270 ♦ a r c dr I .jt + f 



(ia) 



(19) 



(20) 



a r i r 



r 1 ? 

= J 7F (21) 

The pressure rise for the stage is the sum of equations (17) and (19) 

minus the loss due to tip clearance c~ : 

Pt 



C Pf 



Ap f 



r ^ ^ s 



= 2yCU - ^ - a a c dg - Cpt (22) 



r f 3 v f 

A relatively simple expression for the pressure rise in terms 
of the rotor lift coefficient can be obtained by neglecting the 
drag and the tip-clearance terms in equation (22) } which may be 
corrected for afterward, and combining the. result with equation (21) 
to eliminate yo. The result is 
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or 



1 + hi 



c 



f4 + iV-15- - i) 
\W 2 Aa2 C2 2 J 



P f »- 5 J* < xu C H ^ (23) 



where 



c. 



ideal 



B 27aN 



SessSe S risf } (20) «*" n0W be rewitte » *» of the ideal 



7 t 




(24) 



2«r 2 q f 



+ 0„c A 




The expressions (23) and (24) suffice to determine the rotor blade- 
elenent lift coefficient and the relative stream angle at any radial 
location from a desired "frictionless" pressure rise (4p f /q f ) 
Both of these formulas are illustrated graphically in figures fand 5. 
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Figure 4(a) is based on equation (23) and gives o pf against N 
with a r cj as parameter. Figure 5(a) is based on equation (24) 
and gives Cp f . against N with w r /V f or the corresponding 0 r 

a3 parameter. .Figures 4(b) and 5(b) extend the range of the charts 
to high values of N. For convenience , pressure rise is given in 
terms of the related coefficient 



c T 



Ap, 



^ Pf i \pn 2 D 2 

where 



2 

o, 1 - 



p f. - 2 ; 
i 2N 

This relation is plotted in figure 4(b). These charts are useful 
in making quick estimates of fan performance at various design 
operating conditions a3 well as in actual blade design. 

For example, suppose the speed of rotation and axial velocity 

W-y» 

of flow are such that N = — * 1 for a rotor blade element in a 

stator-rotor arrangement. If the blade element is operating at a 
lift coefficient of 0.7 and the annulus is of solidity 0.8, then 
Oci = 0.56 and figure 4(a) gives (Apf/9.f ) ideal = °- 87 - If an 
efficiency of 85 percent is assumed to apply to the pressure, then 
Ap.p 

— =■ = 0.87 x 0.85 = 0.74. Figure 5(a) gives the mean relative 

stream angle 0 r . For N = 1 and (Ap f /q f ) idQal = 0.87, 

0 =39.5°. The blade element should therefore be set sufficiently 

in excess of 39.5° to operate at a lift coefficient of 0.7. 

(See section M Airfoil Characteristics for Cascade Operation.") 



At the design point the ri3e in total pressure (in this case, 
static pressure) is to be constant over the blower plane. Inasmuch 
as the profile-drag terms will not vary much over the blade span, 
the design criterion is taken to be 

c^ = 27GN e ~ - = Constant (25) 
p f,- 9 

11 TtV/ 
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The blade circulation on both stator and rotor will therefore be 
constant along the span. 



The arrangement in which the rotor is upstream of the stator 
will be referred to as a "rotor-stator combination." 

The initial flow into the rotor is assumed to have no rotational 
component, that is, u x » 0 for the rotor. Also, the rotational 
component produced by the rotor is to be completely canceled by the 
stator and, again, 



The only difference between the formulas for this case and for the 
stator-rotor combination is produced by a difference in the sign 
of the stator circulation. Thus, for the rotor, 



Rotor Upstream of Stator 



V 7 * = 7 s a s * 70 



V. 



r 



sin 0 r 




c 




= 27ON - (7a) 2 - 




(27) 



T 



2*r^q f 



(28) 
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r T 



= 2 ? q 

" v r /V f 

For the stator 



<V>J. " J« 



and the pressure rise 



c 



- (ra) 2 - a s c do IS (29) 



Pa w ' 3 Q s V f 

The over-all pressure rise for the stage, including the tip-clearance 
losses c p^j i s 

APf 

= 270S - o r c dr ^ - a s c dB ^ - c pt (30) 

which is the same in form as the corresponding expression for the 
stator-rotor combination. If the friction terms in equation (30) 
are neglected , the frictionless pressure rise may again be expressed 
in terms of the rotor lift coefficient 



18 



NACA AER No. L4Ille 



°r c 7„ " 



It* 




JL pt - ! N 

4N 2 Vi 



2„ 2 



(31) 



and, if equations (26) and (28) are rewritten in terms of the friction- 
lees pressure rise, 



^r 1 
Vf sin ^ 



1 + N 




(32) 



Cm = 



2nr if 



i + Qv^c^ — 



4N j 



Equations (31) and (32) are the basis of the design charts given 
in figures 6 and 7. For the data of the stator-rotor example, that is, 
N = 1 ani OC l = 0.56, figure 6(a) gives (Apf/qf ) ideal = °- 72 and 
figure 7(a) givea a mean relative stream angle of 50.6°. For a fan 
efficiency Of C.85, Ap f /q f = 0.61 in comparison with Ap f /q f = 0.74 
for the stator-rotor arrangement. 
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Rotor Alone 

The initial flow is assumed to have no rotation. The formulas 
for the rotor alone are the same as the rotor formulas (25) to (28) 
pertaining to the rotor-stator combination, except that equation (27) 
for the pressure rise must he rewritten to include the clearance 
losses 



= 2 7 CN - (yo) 2 - CC d i - c Pt (33) 



v f 

or, in terms of the frictionless pressure rise, 




°Pf = c Pfi \~WJ ~ 0Cd T t - c Pt 

Figures 6 and 7 can also be used in this case if the frictionless 
pressure Cp^ of the figures is reduced by the rotational loss 

(7a) 2 = ( c^ /2N N i 2 in addition to the friction terms to get the 



actual pressure rise. The only modification to the rotor-stator 
example is to subtract (^Vf.j 2 ^ 2 = ^ 72 / 2> l 2 * O- 11 from the 
APf/q.f of 0.61. This procedure gives Ap^/q^ = 0.5. 

It is seen that, for a given flow quantity and rotational 
speed, the stator-rotor combination gives the highest over-all 
pressure rise. On the basis of compressibility and maximum-lift 
limitations, however, the rotor-stator combination will yield the 
highest pressure rise. (See section "Illustrative Fan Designs. 11 ) 

It may be remarked that the design criterion for the rotor- 
stator combination is constant circulation 7CN rather than 
constant static pressure over the blower area because a constant 
blade circulation will produce a stable rotational-velocity distri- 
bution in the downstream flow with no radial velocities. As 
mentioned in the section f 'Fan Losses," the rotational velocity 
does not necessarily represent a loss. Even without guide vanes 
as such, if the flow is sent through a radiator, much of the 
rotational energy may be recovered as pressure. 
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The general design criterion for radial-velocity free downstream 
flow for the three combinations discussed is, then, constant blade 
circulation or its equivalent, constant total-pressure increase over 
the blower plane. 

The preceding formulas and the charts for fan design were based 
on the assumption of initially axial flow. If the initial flow has 
rotation, as in operation behind a propeller, the fan design is based 
on the general formulas (8) to (15). For a rotor or stator operating 
in a flow with initial rotation Uj/V fl the rotor-stator design charts 
can still be used to compute Op. in terns of Gc^ and N if the 

abscissa N of the charts (fig. 6) is replaced by N + (See equa- 

v f 

tions (8), (10), and (11).) The ideal pressure rise c pi is regarded 

as 2ya HJ + Y^j. It may be noted also that, for a stator-rotor in a 

flow with initial rotation U]/Vf , the stator-rotor charts can be 
used directly if the design criterion of no downstream rotation 

. ^ f. ¥a • Vr - . .^ r. 

is used. The over-all pressure rise without friction is obtained in 
this case by adding the initial rotational kinetic energy (u-j/Vj.) 2 
to the over-all ideal pressure rise c ?i * 27 r O r N. Although small 

amounts of initial rotation will not greatly modify the fan performance, 
the required blade settings may be appreciably altered. 

The calculation of a multistage fan is carried out on a step- 
by-step basis, each stage being regarded as a single-stage fan 
operating in an initial flow that is the outgoing flow of the 
preceding stage. In order to allow for the successive compression 
of the air by each stage, the fan area may be made successively 
smaller or the successive stages may be redesigned for smaller rates 
of flow. The temperature rise across each stage due to adiabatic 
compression and losses Should also be computed as it may lead to 
excessive operating temperatures in the later stages. 



NACA ARE No. L4Ille 



21 



AIRFOIL CHARACTERISTICS FOR CASCADE OPERATION 

In the preceding sections, the performance of a fan "blade 
element has been related to the lift and drag coefficients at 
which the blade is operating. If the isolated airfoil section 
characteristics (c^ and against a) could he used, the 

knowledge of the blade-element angle of attack relative to the 
mean flow velocity would give the blade- element setting to pro- 
duce the desired performance. This procedure is ordinarily 
followed in the case of propellers for which the solidity is low. 
For high-pressure and hence high-solidity fans, the effects of the 
cascade interference and of the adverse pressure gradient across 
the fan on the blade section characteristics are no longer negli- 
gible. 

The problem of predicting cascade-airfoil characteristics 
for potential flow from a knowledge of the isolated- airfoil 
characteristics has been attacked in two ways. In the first 
method, the isolated airfoil is presumed operating at a certain 
lift and setting relative to the free stream and the necessary 
changes in airfoil shape and setting to produce the same lift when 
the airfoil is operating in a cascade are derived. By the second 
method, the lift of a given airfoil at a given setting in a cascade 
is determined. 

From the first point of view, Betz (reference 5) has given a 
method based on replacing a thin airfoil by vortices in order to 
calculate the flow disturbance of the airfoil. As an example of 
his method, .he takes a 5-percent -maximum- camber circular-arc thin 
airfoil -operating at. c<* =0.7 and places it in a cascade of 
solidity o = 2.03 and mean flow inclined to the cascade at 45°. 
In order bo meixrtain the same lift coefficient of 0.7 in cascade 
operation, the saarber must be increased about 30 percent, with 
relatively sinklJ changes in the position of maximum camber and 
the angle of attack relative to the mean flow. 

. According to the second method, Collar (references 6 and 7) 
has solved for the potential flow through cascades of certain air- 
foil sections including the flat plate.. The conclusion is that 
the relationship from the theory for isolated two-dimensional wings 

q 1 = a Q sin((3 - 0) (34) 

where (3 is the blade-element angle, 0 is the mean-flow angle, 
and (3 - 0 = a is the angle of attack of the blade element, still 
holds provided the right-hand side of equation (34) is multiplied 
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by a factor f and the blade angle 3 with respect to the mean flov 
is reckoned from a slightly different angle of zero lift from that of 
the isolated airfoil. The airfoil shapes given by Collar for which 
the factor f and the cascade angle of zero lift may be derived 
theoretically do not resemble very closely the shapes used in practice. 
A suggested procedure is to use the Collar factor f for flat plates 
in cascade and the angle of zero lift of the isolated airfoil, cor- 
rected for thickness by a method of reference 2. Theoretical curves 
of the factor f as a function of the solidity a = - and blade 

8 

angle p are given in reference 6 and are reporduced in figures 8 
and 9. 

The thickness correction (reference 2) is a change 5 in the 
angle of mean air flov relative to the blade element. It is obtained 
from 



| Da 2 sin (3+0) 
tan 5 = _§ (35) 

1 - | DO 2 cos (p + 0) 



where 

p _ S _ Profile area 
c 2 (Profile chord) 2 

is the reciprocal of the blade profile aspect ratio. The correction 
angle 8 has the effect of increasing 0 of equation (34) to obtain 
the setting for a given lift coefficient; hence, equation (34) now 
reads 

c x = fa Q sin ((3-0-5) (36) 

This equation applies to rotors and to stators downstream of rotors 
if potential flow is assumed. For stators upstream of rotors 
(turbine action), the sign of the right-hand side of equation (36) 
is reversed and angles are measured from the blower plane to the 
convex side of the blade (fig. 10). Note that 0 is the angle 
between the mean flow w and the blower plane. The solution of 
equation (36) for the blade angle (3 involves several successive 
approximations because the factor f and the correction 5 both 
depend on p. 
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Instead of regarding each, airfoil of a cascade as operating at 
a certain lift (and. drag) coefficient with respect to the mean rela- 
tive flow, the cascade as a whole may he viewed as a device for 
turning the entering flow through a certain angle (references 4 
and 8 and unpuhlished work by Quentin Wald of United Aircraft 
Corporation). For example, in reference 4, the investigation of 
a cascade of low-drag airfoils with an initial relative-flow angle 
of 45° and varying exit flow angles, obtained by varying the blade 
angle , led the authors to conclude that r 'the angle through which 
the air is turned by a cascade of blades with a solidity of 1 and 
a small camber is nearly equal to the angle of attack (with respect 
to the entering air) of the blades minus the angle of attack for 
zero lift of the isolated airfoil." 

It may be noted that the turned angle of the flow 0-j_ - 0£ 
is expressible in terms of the blade lift coefficient by means 
of equations (8) and (10), and if drag is neglected, 



2 (cot - cot 0 2 ) 




If the profile drag is large, this expression should be modified 
by the effect of the tangential component of the drag force in 
turning the flow. 

In addition to the mutual interference previously discussed, 
which may be termed a velocity interference, the pressure rise 
(blower action) or drop (turbine action) across the stage affects 
the boundary -layer phenomena in the stage. For example, the 
adverse pressure gradient across a fan rotor reduces the lift 
on the rotor blade, increases the drag, and causes an earlier 
stall. Experiments on a typical airfoil section (reference 3, 
figs. 66 and 67) showed that, for a blade angle (3 of 20° and 
solidity a of 0.85, the factor f of equation (36) was unity 
(although the Collar lift factor f for flat plates was 1.6 in 
this case), 6 was about 1°, and the section stalled at the 
comparatively low lift coefficient of slightly over 0.7. The 
value 0.7 has further been approximately substantiated by 
applying the design: method in reverse to some experimentally 
obtained maximum pressure rises; for example, those of refer- 
ences 4 and 9. 
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On the "basis of the preceding results, which represent typical 

conditions in a high-speed rotor, the factor t 9 the angle 5, and 

the maximum lift coefficient c, may be assumed to he 1.0, 0.0°, 

6 max 

and 0.7, respectively. These values are, of course, provisional and 
subject to change following more extensive experimental investigation. 
For this reason, the rotor blade thus designed should be adjusted in 
operation to give the maximum pressure rise without stalling. If the 
adjustment is small, the error introduced into the blade twist by 
this procedure is probably not serious. The interference corrections 
f and 5 should be included in the design of stators for which 
analyses based on potential flow would be more likely 'to apply. The 
maximum lift coefficient of a stator blade likewise should be sub- 
stantially greater than for a rotor blade. 

As to the' choice of airfoil, an isolated' airfoil section that 
is efficient at the design point Cj may be chosen ( inasmuch as 
there are no comparative data on airfoils for cascade operation) and 
warped by the procedure of reference 5 for cascade operation. 
According to the previously quoted result of reference. 5, a 
15-percent increase in camber of the isolated airfoil section 
should suffice for a cascade of solidity 0=1. If large flow 
deflections are to be accomplished, as by a stator, the choice of 
camber for a stator airfoil section may be guided by the angle 
through which the flow is to be deflected; that is, the camber may 
be determined , by lining up . the tangents at the. leading- and the 
trailing edges of the mean line with the initial and the final flow 
directions, respectively. It is apparunt that more experimental 
results on airfoils in cascade operation are needed to establish 
the. basic data of axial-fan design. 



ILLUSTRATIVE FAN DESIGNS 

The discussions of tho preceding sections will be illustrated 
by the design of a stator-rotor combination, rotor-stator combina- 
tion, and rotor operating alone at an altitude of 40,000 feet. The 
effect of change of altitude will then be considered by estimating 
the performance of the rotor-stator combination at 20,000 feet. 

Each design will be made for maximum pressure rise without 
stalling the rotor bladee and withput running into compressibility • 
effects. The first limitation applies to the hub section because 
the lowest relative velocity, and hence highest lift coefficient 
for a given pressure rise, occurs there. The second limitation 
applies to the tip section because it travels at the highest rela- 
tive speed. These two conditions are expressed as follows: At the 
hub section, by equations (23) and (31), 
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( a r°l r ) 



°*1 



(37) 



max j 2 

At the tip section, by equations (24) and (32) , 





(38) 



R). COR . 

where x^ = — and N Q = 2* The upper sign refers to the stator- 

rotor combination and the lower to the rotor- stat or combination and 
the rotor alone . 

In these equations the values of x h and Vf are fixed by 
the rotor dimensions and the desired flow quantity. The maximum 
value of O r c^ is, or should be, known from experiment, as should 

also w r from the consideration of compressibility. The two equa- 
tions can then be solved for N 0 and c p fi > which determine the 

rotational speed and the ideal pressure rise of the fan. The rota- 
tion parameter N is then fixed for any other location which, 
together with the design criterion of constant c Pf^> determines 

the lift coefficient for any other blade element by using 

either the design charts (figs. 4 and 6) or equations (23) and (31). 
The rotor design is completed by estimating the various losses and 
finally the blade angles to give the required lift coefficients. 
The stator, if any, is designed from a knowledge of the ideal pres- 
sure rise and the condition of complete cancelation of stream rota- 
tion. The details are given in the following discussion. 

The solution of equations (37) and (38) for c p{ ^ and N Q 
is obtained by solving the quadratic 



au 2 + bu - c = 0 
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as 



u = 



+ 4ac 



2a 



(39) 



and then substituting 



N 0 - a + u I 

c p = 4N c u 
1 J 



(40) 



where 



> 



m = .1 



a = 



Xh(OrCi r j, 



- 1 



max 



c = 



1 - " *T 

2mk 2 6- - x h 2 ^ 
k 2 (l + % 2 m 2 ) 



> 



(41) 



The upper sign refers to the stator-rotor combination and the lover , 
to the rotor-stator combination and the rotor alone. 

As to the fan dimensions, the outer radius will probably be 
fixed by considerations of engine or cowling diameter. Assume 
R 0 m 2 feet. The inner radius should be such as to keep the solidity 
at the hub down to a reasonable value and not too different from the 
solidity at the tip. Also, the resulting fan area should be large 
compared with the tip-clearance area and should permit efficient 
diffusing downstream of the fan. Assume an inner radius of 1.55 feet 
giving a radius ratio x h = 0.776, a fan area Af = 5 square feet, 
and a blade span of 5.375 inches. If a tip clearance of l/8 inch 
or 2.3 percent of the blade span is further assumed, then, according 
to figure 3, a loss of 8.3 percent of the ideal pressure rise will 
occur. 



NACA ARE No. L4Ille 



27 



The maximum lift coefficient is taken as c 7 =0.7 in view 

6 max 

of results quoted previously from reference 3. According to some 
results of reference 2, a rotor with "blade chord increasing toward 
the outer section is more efficient than one with either constant 
chord or chord decreasing toward the outer section. The reason is 
connected with the fact that the secondary flows previously 
mentioned (see section "Theory of the Axial-Flow Fan 1 ') cause a 
dead-air flow toward the blade tips, which may consequently stall 
first unless the lift coefficient is kept. low. A constant- solidity 
blade with <j r = 0.85 is accordingly chosen. The quantity ayCj 

is therefore 0.7 x 0.85 = 0.595. 

The maximum relative velocity is determined by the somewhat 
arbitrary though conservative compressibility criterion that the 
maximum local relative velocity at the rotor blade is not to exceed 
0.9a, where a is the speed of sound at the fan. Inasmuch as the 
maximum local relative velocity is about 1.31 times the relative 
velocity and the speed of sound at 40,000 feet is 1000 feet per 
second, the maximum relative velocity is 

Wr max = i%| x 1000 = 686 feet per Qecond - In order to get V f , 
the required mass flow is taken as M = 0.821 slug per second 
(see section "Power Losses in a Cooling Installation") which, with 
a density at the fan of 0.000616 slug per cubic foot and the fan 
area of 5 square feet, give3 Vf ■ 266 feet per second. Therefore, 
(*rAf) max = 2.58. 

The values of x h , (a r c? ) , and (w r /V> ) thus deter - 
n max r 1 'max 

mined are substituted in equations (41), (39), and (40) to give 

the values of N D and Cp f ^ shown In table I. The values of 

rotational speed derived from HT 0 are also given. These values 
of N 0 and c pf # are the basis for the various fan designs given 

in tables II, III, and IV. The tables contain the itemized design 
procedure and the necessary references to the text and figures. 
Some of the results aro included for comparison in table I. It 
appears that the rotor- stator arrangement is capable, of producing 
the greatest pressure rise, though \cith the requirement of operating 
at a considerably higher rotational speed. From a consideration 
of the solution of equations (37) and (38), this result is found 
to be generally true if 



m2 (l - 2x h 2) < 1 + [HWP p . % 2) 2 ( m 2 . y 
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Because the second term on the right-hand side is quite small and 
usually positive, the rotor-stator combination is seen to yield the 
highest pressure rise if x h = Rh/K 0 £ 0.707. BeJLow this value the 
rotor-stator or the stator-rotor combination can yield the higher 
pressure rise, depending mainly on the value of m. 

For operation at an altitude of 20,000 feet, the essential dif- 
ference is in the air quantity required for cooli.ng. It is assumed 
that, at 20,000 feet, 1.25 times the mass flow anid hence 0.475 times 
the volume rate of flow at 40,000 feet are necesetary. If the rota- 
tional speeds wore maintained at the original val.ues, the rotor 
blades would be completely stalled. As can be se.on from the perform- 
ance curves in references 2, 9, and 10, stalling the blades means 
that the maarlmtlB pressure rise is developed but with high losses. 
If the rotational speeds :rre diminished in the ratio of 0.475, then 
N is the sane and the nondimensional operation o f the fan is the 
same as at 40,000 feet. Therefore, because 



V 2 /20,000 V 2 / 



40,00(3 

the actual pressure rises of any of the combinations will be about 
0.6 of the values at 40,000 feet. With fixed blades this performance 
is about all that can be achieved at the lower altlitude on the 
assumption that c ly , =0.7. If the rotor and s tat or blade angles 

are adjustable, however, and the rotor speed is maintained at 3310 rpm 
for the rotor-stator arrangement, then for a hub-section value of 

O r c 7 = 0.595 and N h = |4H " 4 ' 26 > figUre 6(b> giVGS * 
r 4 r 0.4/5 , 

(Ap f /q f ). of 9.8. If the same efficiency as art 40,000 feet is 

assumed, the actual Ap f /q f is 0.9 X 9.8 = 8.8 and Ap f is 
22 inches of water, as compared with 9.1 inches at 4:0,000 feet. The 
twist of the blade is incorrect at the lower altitude, however. In 
this case it is desired to calculate the pressure rirae at any radial 
location from a knowledge of the blade angle and N for the location. 
The solidity a and interference corrections f and 6, if any, 
are also known. The ideal pressure rise is given directly m terms 
of these quantities by eliminating c Zr from 
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and , by equation (36), 



Cj = fa 0 sin (6 - 0 r - 8) 



where 




The result is 



c 



OrfaQNptf sin (p - 8) - cos (p - 5)] 
p « _ 

1 ± JL_£ sin (p - 5) 
4 



(42) 



where the minus 3ign refers to the stator-rotor combination and 
the plus sign refers to the rotor- stator combination, which includes 
also the case of rotor alone. (See section "Kotor Alone.") It is 
assumed in equation (42) that, if a stator is used, it is kept 
adjusted in such a way that no stream rotation exists downstream 
of the "blower. The error' caused by using equation (42) with fixed 
stator is small for a rotor- stator but not for a stator-rotor 
combination. In the stator-rotor arrangement, the stator and rotor 
should be analyzed separately by means of the general formulas (8), 
(10), (11), and (36). The actual pressure rise is obtained from 
equation (10) applied to rotor and stator. This pressure rise will 
not actually be obtained if the rotor maximum section lift coef- 
ficient is exceeded. Furthermore, if the pressure rises vary 
greatly over the blower area, they should be averaged by the method 
of reference 2. 



The design method presented in the preceding sections has been 
used to compute the known performance of existing fans such as that 
of reference 9. The results of these and of other computations for 
new designs, which have subsequently been tested, show that the 
method predicts performance satisfactorily, although an experimental 
adjustment of the blade setting may be necessary. (See section 
"Airfoil Characteristics for Cascade Operation. 1 ') 
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POWER LOSSES IN A COOLING INSTALLATION 

It is appropriate to conclude vith a discussion of the internal- 
cooling-power losses in a cooling installation that includes a fan. 
The statement in the "Introduction 11 that the pover required to drive 
a fan may be more than recovered by means of the higher exit velocity 
of the cooling air thus attainable will incidentally be justified. 

The internal cooling power is the sum of the power input to the 
fan and the momentum power of the cooling air. The power input to 
the fan per unit mass flow of cooling air is a function of the fan 
pressure rise and the state of the air in front of the fan. The 
momentum power per unit mass is the product of the free- stream 
velocity and the difference between free- stream and ultimate wake 
velocity of the cooling air. External cooling power , such as caused 
by excessive flap deflections with resulting increase of drag, is 
not discussed herein. This component of the total cooling power 
can be included in an analysis of an actual cooling installation 
by means of experimental cowl-flap data. Likewise, drag increments 
due to the weight of the component parts are not considered. 

A cooling installation of cowling, fan, cooling resistance 
(such as engine or radiator), and exit is represented schematically 
in figure 11(a). Cooling air ; regarded as compressible, enters and 
undergoes adiabatic but not isentropic compression from state 0 
to state fl. , The total energy per unit mass of air remains 
constant. in this process, that is, 

Vo 2 Vfi 2 
c P T o +■ ~ 1 c P T fl + ~ 

but the inlet duct causes a stagnation pressure drop Ap^. By a 
stagnation pressure drop across a resistance is meant a pressure 
drop measured between, or corrected to, two stagnation regions 
before and behind the resistance. Next, the fan increases the 
total energy per unit mass by 

§=c p (T f2 -T fl ) + |(V f2 2-V fl 2) . (43) 

while the pressure rises from p fl to p f2 . The air then passes 
through the fan diffuser in which a stagnation pressure drop Ap d 
occurs at constant total energy. At state 1, which is the near- 
stagnation condition just in front of the cooling resistance, a 
fraction of the engine-heat rejection may be added to the part of 
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the flow that goes through the resistance proper. This flow then 
expands isentropically into the resistance area, state 2. From 
state 2 to the resistance exit, state 3e, the cooling air is 
increased in total energy by the remainder of the heat rejection 
and suffers most of the cooling pressure drop. At state 5e the 
flow is ejected into the large space at state 3, the pressure and 
total energy remaining constant, with the result that the kinetic 
energy of the air at state 3e is largely lost. The cooling- 
resistance pressure drop is thus Ap r = P3 - Pj. The entire flow 
at state 3 then expands almost isentropically, with constant total 
energy, from its practically stagnation condition at state 3 to 
the ultimate exit conditions p 0 ' and V 4 . 

The power input to the fan is given by equation (43). The 
second term on the right is negligible for a single-stage fan with 
constant flow area, whereas in a multistage fan the flow area is 
adjusted to keep a constant axial velocity. The first term can be 
written in terms of the initial state and pressure ratio across 

Zli 

the fan by using the isentropic compression formula T c: p 7 . 
The power input to the fan thus becomes 

P f = Mc p (T f2 - T fl ) 



7MP! 



(7-i)n f p 1 




MAp f 

T) f Pl 



1 - _L . Jl£ + 
27 v x 



(44) 



The fan efficiency r\f in this expression takes account of the 
fact that the compression is not quite isentropic and that conse- 
quently the actual final temperature Tfo is greater than the 
isentropic final temperature would be for the same pressure rise. 

The momentum power of the installation is, by the momentum 
equation, 



(45) 
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The ultimate exit velocity V 4 , at vhich the pressure is that of 
the free stream p Q , is determined by the condition of constant 
total energy from states 3 to 4: 



c p T 3 ■ c p T 4 + 



T 4 




(46) 



The quantity 1 - -3 is the Meredith effect efficiency, that is, 
T 3 

the efficiency of conversion of energy at state 3 to kinetic energy 
of the wake at state 4. For isentropic expansion between states 3 
and 4, it may be expressed as 



T 3 



P3 



(47) 



If conditions between states 3 and 4 differ considerably from the 
isentropic, equation (47) must be modified. For well-designed out- 
let ducts, no modification of equation (47) is necessary. The pres- 
sure at state 3 is the free-stream stagnation pressure p a plus the 
sum of the stagnation- pressure changes up to state 3 



?3 ■ Ps + A Pf " Ap i " Ap d " Ap r 



(48) 



The total energy per unit mass at state 3 is the free- stream total 
energy plus the total-energy additions in the installation: 



V 0 * Pf H 
c p*3 = c P T o * — * M" + M 



(49) 



The substitution of equations (49), (48), and (47) in equation (46) 
yields for the exit velocity 



2 - 1 P ° + 2 + M + M 



Po 



\ 7 



v p- + &V* - Ap, - Ap, - Ap 



(50) 
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which determines the momentum power by equation (45). The sum of 
equations (44) and (45) is then the internal-cooling-power loss 
of the installation. 

Before equation (50) can "be used, the cooling and fan-diffuser 
stagnation- pressure drops Ap^ and Ap^ as well as the required 
mass flow M must be determined as a function of the fan pressure 
rise Ap~. This procedure is necessary because the fan ahead of 
the cooling resistance increases the temperature and density 
ahead of the cooler, thereby increasing the required mass flow 
for cooling but decreasing (not always) the cooling pressure 
drop. The determination is made by a step-by- step process similar 
to the processes of references 11 and 12. 

In order to get M, Ap r , and Ap^ as functions of APf, 
the pressure and temperature just ahead of the cooler are needed. 
The following equations determine the pressure and the tempera- 
ture in front of the cooler, as well as the mass flow, as func- 
tions of fan pressure rise: 
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(a) Free stream to state fl, constant total energy: 



c p T G + |V Q 2 = c p T fl + |V fl 2 



(b) Pressure and density at state fl, inlet diffusing 
efficiency tj.j 



Pfl 



- (i - niho 



Pfi » V RT fi 



(c) Pressure and temperature rise across fan, fan 
efficiency r\f. 

__2L 



Pf2 
P fl 



T 



f2 



n f — + (i - n £ ) 



7-1 



(d) Pressure and density at state f2: 

/Pf2\ 



P f2 ■ Pfl| 



p fi; 



A Pf . p f2 



p f2 = P f2/ RT f2 

(e) State f2 to 3tate 1, constant total energy: 

c p T l ' = C P T f 2 + fyt 2 

(f) Pressure at state 1, fan-diffuser efficiency X\^i 



V, = P 



f2 



where 



_2L 
V 7=I 



- (i - 



d'^f2 



1f2 = |pf2 V f 2 2 



(g) Variation of mass flow M with cooler entrance temperature 
Tj_', where Tj_' is the temperature in state 1, before any 
of the heat is added to it. This variation is determined by a 
cooling correlation for the cooler. (See references 11 and 12.) 
Figure 12 illustrates this variation for the air-cooled engine 
of the illustrative example given hereinafter. 



(h) Equation of continuity at fan 
M = Pf2 v f A f2 



M = PfiVfA fl 
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In these equations the atmospheric conditions p Q; T Q , flight 
speed .V 0 , and efficiencies r\^ f ; are all known. Further- 

more } either the fan entrance area Afj_ is given and the fan exit 
area Afg regarded as variable to maintain constant axial velocity 
Vf, or vice versa. In the following illustrative example, tile 
value of A-q was assumed ; although about the same results would 
bo obtained in either case. With an assumed Vfy the following 
quantities are determined in order from equation (5l): Tfj_ from 
(a); Pfn and p fl from (b); M from (h); from (gj; T f2 

from (ej; P f2 /Pfi from (c); p f2 , Ap f , p f2 from (d)j q f2 
and p-j. from (f). The temperature is then increased to 

by the heat rejection, if any, ahead of the cooler. In addition, 
the value of G = p£V 2 is determined as a certain fraction of M. 
The proportionality between G and M can be obtained from the 
sea-level data for the cooler by performing the cooler pressure- 
drop calculation in reverse to find G corresponding to the 
known pressure drop. 

From the values of pj_, Ti, and G corresponding to a fan 
pressure rise Ap^ , the pressure drop Ap r across the cooler can 
be determined in two steps. First, if only the air that goes 
through the cooler proper (in an air-cooled engine this amount is 
less than the total mass flow) is considered , the statu of the air 
at station 2 is determined by Bernoulli's equation, which may be 
written 



pi 



7 



p ? 1 1 R Z C, 2 T-> 

- S K (52) 



7+1 Boj« 



g) 



7 



This equation can bo solved for p 1( /p 2 with the aid of figure 13, 
which gives Pj/p 2 against K for 7 = 1.4, B = 1715, and 
Cp = 6020 in foot-pound-second units. Thu temperature at 
station 2 is' given' by 

Hi 



T i fnV 

p 2 , 



T 2 
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Secondly, the pressure at the baffle exit, station 3e, is determined 
by an equation based on heat-transfer and skin- friction considerations. 
(See references 11 and 12.) 



log Q -£ 



P2 T v 



*2 ?3e 



P2 £ 



P3 T w "T 3 e T 2 ~ RG 2 (T 2 + T 3e ) 



P2/ 



(53) 



in which T 3e is greater than T 2 by an amount corresponding to 

the heat added from stations 2 to 3e and T v is the estimated average 

engine temperature. Equation (53) is plotted in figure 14 for a range 



of the parameters 



T. 



v/ 



T 



3e 



and 



P2 



T v - T 3e T 2 



EG 2 (T 2 + T 3e ) 



The total energy c p T 3 behind the cooler is calculated by equa- 
tion (49), using the first one of equations (44) for the Pf/M term. 
The values of c p T 3 and p 3 can now bo substituted directly in 



equation (46) to give the exit velocity 
completed by equations (44) and (45). 



and the calculation is 



For the fan behind the resistance and for the station numbering 
as in figure 11(b), the expression through equation (50) previously 
derived hold without change. The calculations are greatly simplified 
in this case because M and Ap r are constant at the values given 
by the preceding method for Ap f = 0. Conditions at the fan can be 
conveniently determined from the stagnation conditions behind the 
cooler and the mass flow by using figure 13. 



The calculations previously outlined for the fan ahead of the 
cooling resistance were carried through for an air-cooled engine. 
The following data, based on an actual proposed cooling installa- 
tion, were assumed: 



Altitude (Army air), feet . . 40,000 

Atmospheric pressure, p 0 , pounds per square foot 391.9 

Atmospheric density, p Q , slug per cubic foot 0.000550 

Atmospheric temperature , T Q , °F absolute 415.4 

Specific heat at constant pressure, c p , foot- 
pounds per slug per °F absolute 6020 

Ratio of specific heats, y 1.40 

Inlet-duct efficiency, y\ ± 0.90 

Fan efficiency, T) f • 0.80 

Fan-diffuser efficiency, r) d 0.90 

Fan- entrance area, &fi> square feet 5 

Engine heat rejection, H, foot-pounds per second . . . 550,000 
Flight speed, V Q , feet per second 587 
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In addition the following assumptions were made: Figure 12 for M 
against f ; S = 0.305M; .one-half the heat rejection added ahead 
of the "baffles and one-half in the baffles; 80 percent of the total 
mass flow through the "baffles; the average engine temperature 
T v a 8880 y absolute. 

The results are given in figures 15 and 16. Some of the 
various cooling- air quantities are shown in figure 15. The varia- 
tion with fan pressure rise is practically linear in all cases , a 
fact that may be used to reduce the labor of subsequent calcula- 
tion. The engine pressure drop p 3 - v± decreased from 26 inches 
of water at Ap f = 0 to 24 inches of water at Ap f m 38 inches of 
water. For high-enough fan compression, the temperature of compres- 
sion tends to compensate 'for the effect of the increased density in 
reducing the engine pressure drop. 

The variation of fan and momentum power and of their sum, the 
internal cooling power, is shown in figure 16. The internal cooling 
power is large, varying from a maximum of 640 horsepower to a mini- 
mum of 330 horsepower. It is seen that, without a fan, there is a 
cooling-pressure deficiency of 8 inches of water. The cowl-flap 
deflection required to overcome this deficiency (about 0.5q o ) 
would result in a cooling-power loss much greater than the indi- 
cated maximum value of 640 horsepower. A fan is therefore 
essential. The maximum estimated pressure rise from a single- 
stage fan is indicated in figure 16 as 8.5 inches of water. This 
value was obtained by multiplying the value of 9.1 inches in 
table I by . 80/85.6 to make the fan efficiency consistent with 
that assumed in this example. This pressure rise reduces the 
internal cooling power to 520 horsepower. Two stages would decrease 
this value to 350 horsepower, which is only 20 horsepower above 
the minimum possible cooling-power loss. In order to reduce the 
minimum cooling-power loss, the required engine pressure drop for 
the given heat rejection would have to be reduced, for example, 
by improving the finning. The flight speed chosen for the calcu- 
lations was 400 miles per hour. If the speed is reduced, as in 
climbing, the cooling-pressure deficiency increases. For example, 
at 310 miles per hour, the deficiency would be about twice the 
value at 400 miles per hour. A two-3tage fan would therefore be 
necessary to enable climbing without opening the flaps. 

It is of interest to note that, at the point of minimum total 
internal power, the fan power is increasing at the same rate as 
the thrust (negative momentum) power; hence, the efficiency of 
conversion of additional power input to the fan into thrust power 
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is 100 percent at this point. This efficiency drops off for larger 
additional power inputs to the fan. The possible saving of power , 
as compared with putting this power into the external propeller, 
however, is small. 

The question of the advisability of locating the fan behind 
the engine, or cooling resistance, can be answered in the negative, 
at least for conditions approximating those of the example. Some 
comparative calculations of fan before and behind the engine showed 
that the cooling power was of the order of twice as much as with the 
fan in front. Furthermore, the maximum performance of a single- 
stage fan was only about 60 percent of that with the fan in front. 
The main reason for this behavior was that the density behind the 
engine was about half that in front. Consequently, the volume rate 
of flow for a given mass flow was double and the power input to the 
fan for a given pressure rise, roughly QAp, about double. A con- 
tributing factor is. the increased engine pressure drop, which shifts 
the momentum-power curve to the right and consequently increases 
the total internal power. The decreased maximum performance of a 
single stage is a consequence of the lower density in combination 
with the higher ratio of axial-flow velocity to fan rotational 
speed. 

The method of internal-power calculation given is an accurate 
one. The results are sensitive, however, to the values of the 
various pressure drops, particularly to the pressure drop of the 
cooling resistance. It is important, therefore, to be reasonably 
sure of the data that determine the engine pressure drop before 
indulging in a detailed internal-power analysis. It may be 
mentioned here that the method of calculating the exit velocity 
by Bernoulli f s equation for incompressible flow, that is, 



can lead to large errors on the pessimistic side, particularly for 
large pressure and density changes through the system. 

It is evident that the choice of an axial-flow fan for high- 
altitude cooling depends on considerations of internal cooling 
power as well as on marginal cooling-pressure deficiency, struc- 
tural, and mechanical considerations. In the illustrative example, 
a two- stage rotor- stator fan arrangement ahead of the engine would 
be recommended. If the rotational speed could not be achieved 
because of mechanical difficulties, a stator-rotor arrangement 
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could be used although at the probable sacrifice of some performance, 
(See table I.) The stator-rotor arrangement might also be advan- 
tageous if the flow conditions ahead of the fan were not uniform, 
with the result that the guiding and stabilizing effect of a set 
of stator vanes was needed. Some type of rotational- speed adjustment 
is also desirable , particularly on high-speed fans, in order to avoid 
blade stalling at the lower altitudes where smaller volumes of air 
are required. 



CONCLUSIONS 

The following conclusions are indicated by the results of the 
present study of the theory, design, and effect on cooling of the 
constant-velocity axial-flow fan as related to the problem of 
high-altitude aircraft engine cooling; 

1. An axial-flow fan may be used to give a moderate pressure 
rise to a large flow quantity of air. Typical performance of a 
4-foot-diameter single-stage rotor- stator fan at 40,000-foot alti- 
tude is 9 inches of water pressure rise at a flow quantity of 
80,000 cubic feet per minute and a rotational speed of 3300 rpm. 

2. Existing theory provides a satisfactory basis for the design 
of axial-flow fans. 

3. The use of stationary guide vanes and small rotor tip 
clearances is essential for efficient fan design. 

4. Further research is needed for the development of airfoil 
sections especially suitable for use in high-speed, high-solidity 
rotors. 

5. The losses in internal cooling power at high altitude are 
large when the cooling is marginal. 

6. An axial-flow fan more than pays for itself in cooling 
power saved when cooling conditions are marginal. 

Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I - FAN PERFORMANCE BASED ON 
MAXIMUM-LIFT AND COMPRESSIBILITY 
CRITERIONS 



Fan 
arrange- 
ment 


coR Q 
N_ = —2 

° V f 


+ 


Rota- 
tional 
speed 
(rpm) 


Pressure 

rise 
(in. of 
water ) 


Input 
(hp) 


Effi- 
ciency 
(per- 
cent) 


Stator- 
rotor 


2.12 


2.17 




2700 


8.1 


120 


85.5 


Rotor- 
stator 


2.61 


2.40 


3310 


9.1 


133 


85.6 


Rotor 
alone 


2.61 


2.40 


3310 


7.9 


133 


75.1 
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TABLE II - ILLUSTRATIVE FAN DESIGN - STATOR -ROTOR COMBINATION 
[NACA 6512 blade section] 



Fan quantity 


Source 


Calculated 






values 




Rotor, 32 blades 






0.776 


0.888 


1 


r, ft 




1.552 


1.776 


2 


i 




1.649 


1.885 


2.122 




Fig. 4 


0.595 


0.481 


0.396 






0,85 


0.85 


0.85 












G r and 32 blades 




3.11 


3.56 


4.01 

• 


cj r for constant Cp^.. 


/Fig. 4 or 
v equation 


0.7 


0.566 


0.466 


v r /V f 


tM f-* cz 

Fig. 5 




9 


C • DO 


\ * / 4 J- 


Equation (40) 


2.170 


2.170 


2.170 




Soe stator 
data 


0.0174 


0.0151 


0.0133 


( c d 0 = 0,013, assumed) 










Equation (19) 


0.0245 


0.0264 


0.0285 


(cd r = 0.013, assumed) 










/Ap f \ 

— i. j for 2. 78-percent 


Fig. 3 


0.180 


0.180 


0.180 


\ ^f /x. clearance 
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TABLE II - ILLUSTRATIVE FAN DESIGN - Continued. 



Fan quantity 


Source 


Calculated 
values 


Rotor, 32 blades 


( 


r APf\ / A Pf"\ 












r^f ) ~ \~^f ). 

K 1 'resultant x A 'i 




1. 948 


1.948 


1.948 




\ q f 4 a A *f / dr A «f/ t _ 












Ap f , in. vater 




8.14 


8.14 


8.14 




Equation (12) 


1.365 


1.207 


1.090 


Section efficiency, 
percent, 












A ^ f /resultant 
Nc T 


Equation (15) 


86.6 


85.6 


84.2 


Output, QApf, hp 






102.5 




Input . 0ut P ut , hp 
Av. efficiency 






120 




Mean flow, deg, 












0 - sin-l — 1 

V r /V f 


Fig. 5 


26.77 


24.71 


22.82 


P 


- 0 with no inter- 
ference corrections f 
and 5 and a Q = 5.65 


Equation (34) 


7.11 


5.74 


4.74 


3 measured from zero-lift 

line 




33.88 


30.45 


27.56 


B 


measured from chord line 
for a 7 =6.2° 




27.68 


24.25 


21.36 
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Fan quantity- 


Source 


Calculated 
values 


Stater, 39 blades 


x = r/R Q 






0.776 


0.888 


1 


O f r.r 

d 4N 




Equation (25) 


0.329 


0.288 


0.2555 


C g assumed to give 
reasonable c 7 






1.275 


1.116 


0.990 


Chord 9 c S; in. from 
and 39 blades 






3.83 


3.83 


3.83 


?s 






0.516 


0.516 


0.516 






Equation (16) 


1.051 


1.0405 


1.031 


C7 = ; 

3 V V f 




T?n i IP *i on f ft I 




u • y ji 


1 
-L 






Equation (17) 


0.0174 


0.0151 


0.0133 


0 = sin" 1 deg 






72.0 


73.75 


75.70 


5. Da 2 for D = 0.08 
6 






0.0682 


0.0524 


0.0411 


f 1 

5 > for a 0 = 5.65 
0 + 5 - B J 




0.58 
2.88 
17.41 


0.62 
2.11 
16.44 


0.66 
1.55 
15.55 


P 






57 .47 


59.42 


61.70 


B' for a 2 = 6.2° 

0 






63.67 


65.62 


67.90 j 
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TABLE III - ILLUSTRATIVE FAN DESIGN - ROTOR -STATOR 

COMBINATION 
[NACA 6512 blade section) 



Fan quantity 


Source 


Calculated 
values 


Rotor 3 32 "blades 






0.776 


0.888 


1 


r, ft 




1.552 


1.776 


2 






2.022 


2.316 


2.609 


a r c Z r 


r Xg. 0 






U . ODD 


a r 




0.85 


0.85 


0.85 


Chords c r , in* from O r 
and 32 blades 




3.11 


3.56 


4.01 


°l r for constant 


/Fig. 6 or 
\equation (31) 


0.7 


0.536 


0.421 


v r /V f 


Fig. 7 


1.994 


2.287 


2.58 


PS), ■ > 


Equation (40) 


2.404 


2.404 


2.404 


(c(i a = 0.013, assumed) 


See stator 
data 


0.0157 


0.0136 


0.0120 


V «f/ d * dr v f 

u r 

(c^ ta 0.013, assumed) 

r i 


Equation (27) 


0.0220 


0.0253 


0.0285 
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TABLE III - ILLUSTRATIVE FAN DESIGN - Continued. 



Fan quantity 


Source 


Calculated 
values 


Rotor ; 32 blades 


/a x 

(— ) for 2.78-percent 
V If clearance 


1 


0.200 


0.200 


0.200 






2.166 


2.166 


2.163 


v x ^resultant 










Apj, in. vater 




9.06 


9.06 


9.05 


Crji 


Equation (28) 


1.227 


1.090 


0.990 


Section efficiency, percent, 










V /re suit ant 


"CVi 1 1 o "1* T r> t» [1m 
JZa^Uo. b XU.ll \XU/ 


P>1 ? 

KJ ( % Lt 


as a 


83.8 


T] = 

' Nc T 










Output, QApf, hp 






115 




Input = Output/Av. effi- 
ciency, hp 






133 




* = 8in " 1 557*? des 


Fig. 7 


30.13 


25.94 


22.81 


P - 0 for a 0 =5.65 


Equation (34) 


7.11 


5.44 


4.28 


13 




37.24 


31.38 


27.09 


6' for * 6.2° 




31.04 


25.18 


20.89 



National Advisory Committee 
for Aeronautics 



NACA ARR No. L4Ille 



47 



TABLE III - ILLUSTRATIVE FAN DESIGN - Concluded. 



Fan quantity 


Source 


Calculated 
values 


St at or, 39 blades 


X 






0.776 


0.888 


1 


7C (Ap f /q f ) i 
2 4N 




From rotor 
calculations 


0.297 


0.260 


0.231 


Q 3 (assumed) 






1.159 


1.015 


0.901 


Chords c s , In; from 
and 39 "blades 


°s 




3.49 


3.49 


3.49 


?s 






0. 513 


0.513 


0. 515 








1.041 


1.032 


1.026 


2-y 

c 

2 B ~ V 3 /V f 




Equation (8) 


0.986 


0.995 


1 


U;j d - CsCd sv f 

a s 




"PYm^t ion f?Q^ 




0 0136 


0 0120 








73.69 


75.51 


76.90 


f > for D s= 0.08 
P-0-8 J 




Fig. 8 and 
equations (35) 
and (36) 


0.0562 

0.525 
0.67 
19.42 


0.0432 

0.58 
0.45 
17.68 


0.0340 

0.63 
0.31 
16.31 


0 






93.78 


93.64 


93.52 


(J ! for a 7 = -6.2° 






87.58 


87.44 


87.32 
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TABLE IV - ILLUSTRATIVE FAN DESIGN - ROTOR ALONE 
[NACA 6512 blade section; 32 blades] 



Fan quantity 



r 
N 



Source 



Ap f \ 

Cp f for rotor- 



> (a) 



W 4 



stator 



If/ 



rotational energy 

Uvf/<if){ 2 



= (7a) 2 = 



211 



Calculated 
value 3 



(a) 



Equation (33) 



(a) 



(a) 



2.404 12.404 



0.353 



0.270 



2.404 



0.212 



a Same as in table III. 
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FAN 'DESIGN - Concluded. 
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Fan quantity 



If/ 

resultant 

Ap^, in. water 

Ctji 

Section efficiency } tj 

Output, Q(Ap f ) 

1 av hp 

Output 



Input = 

0-0 

13' 



Av. efficiency ' 



, hp 



00 



1 



>(a) 



Calculated 
values 



0.022 0.025 



0.200 

1.829 

7.64 
1.227 
73.6 



0.200 

1.909 

7.98 
1.090 
75.6 



0.029 
0.200 

1.963 

8.20 
0.990 
76.0 



100 
133 



(a) 



(a) 



(a) 



^ame as in table III. 

^Formulas and references same as in table III. 
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Fig. 1 
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Figure 2. - Rotor forces and velocities* 




Figure 3 # - Tip-clearance losses in percent of ideal pressure rise. 
(Average of curves taken from reference 2 # ) 




Figure 4 # - Stator-rotor design chart relating ideal 

rise, N, and <5 r 0 l r * 



pressure 
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Fig. 4b 




(b) High N. 
Figure 4.- Concluded. 




(a) Low N. 

Figure 5.- Stator-rotor combination design chart relating ideal 

pressure rise, N, and w r /Vf. 
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Fig. 5b 





(a) Low N. 

Figure 6.- Rotor-stator combinatioa and rotor alone design chart 
relating ideal pressure rise, N, and o^c, • 
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lb) High N. 
Figure 6.- Concluded. 




(a) Low N. 



Figure 7.- Rotor-stator and rot ~>r-alone design chart relating 
ideal pressure rise, N, and w r /V^. 
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Fig. 7b 




Fig. 10 
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Fig. 11a, b 
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(a) Fan ahead of cooler. 




(b) Fan behind cooler. 
Figure 11.- Schematic cooling installation. 




Figure 12.- Assumed mass flow M against air-cooled engine entrance temperature T^ 1 



*3 




(a) Large values of pressure ratio. 



Figure 14. - Graphical solution of equation (53). 



t 




Figure 15.- Coollng-alr condition against fan pressure rise. 




Figure 16,- Internal-cooling-power losses 



against far pressure rise. 



